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The perturbative QCD approach to quarkonium decay into a photon and hadrons is reconsidered. It is shown 
that a consistent treatment within perturbative QCD calls for the introduction of a fragmentation contribution 
which has been neglected so far. The ensuing phenomenological implications are discussed, and, in particular, 
the possibility of measuring the gluon fragmentation function of the photon is addressed. 



1. Introduction 

One of the earliest applications of perturbative 
QCD was the calculation of the decay widths of 
heavy-quark-antiquark bound states (quarkonia) 
[1]. The use of perturbative QCD in this con- 
text is based on the fact that these decay pro- 
cesses involve two very different physical scales: 
the binding energy e of the bound state and the 
heavy-quark mass mq . Since mq ^> e for the \?- 
and T-system, their decay width T into a given 
final state / can be factorized as follows 

T(X Q q^f) = Rl\M Q q^ f \ 2 + 0(v 2 /c 2 ) . (1) 

Here Ro is the wave function of the bound state 
at the origin, whilst M is the matrix element 
describing the annihilation of the on-shell heavy 
quark and antiquark into the final state /. 

As long as the relativistic corrections 0(v 2 /c 2 ) 
in (1) are small, all of the bound-state effects 
are factorized into the wave function, and can- 
cel if one considers ratios of decay widths for dif- 
ferent final states. In this case one is thus left 
with the evaluation of the matrix element M , 
which is dominated by large momentum scales 
(of the order of mq). This term is computable 
(modulo hadronization corrections) as a pertur- 
bative series in as(mq) provided it is infrared 
and collinear safe, that is, the decay process is 
fully inclusive with respect to the final-state par- 
tons. 



A well-known example of decay rate which is 
computable perturbatively is the ratio between 
the hadronic and leptonic widths. It has been 
evaluated in leading (LO) and next-to-leading 
(NLO) order in QCD perturbation theory [2,3] 
and is mainly used for a^-determinations [4-6]. 

In this paper we concentrate on the inclusive 
radiative decay 

quarkonium — > 7 + hadrons , (2) 

where 7 denotes a prompt photon, i.e. a final- 
state photon not coming from hadronic (essen- 
tially, 7r and 77) decays. According to common 
wisdom, the decay rate for the process in Eq. (2) 
is computable in QCD perturbation theory and 
even a NLO result is available in the literature 
[7,8]. Despite common wisdom, we would like 
to point out that this perturbative-QCD pic- 
ture is not correct, at least from the theoretical 
viewpoint. In fact, it (surprisingly) misses out 
a leading-twist and leading-order fragmentation 
component 1 . 

In the following we clarify this issue and con- 
sider its phenomenological implications. A more 
detailed discussion will appear elsewhere [10]. 

2. Xqq — > 7 + hadrons : revised theory 

Let us compare the hadronic decay Xqq — » 
hadrons with the radiative decay in Eq. (2). In 
the first case the corresponding partonic final 
state consists of all the possible channels with 
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1 As far as we know, the need to include a fragmentation 
component was first pointed out in Ref. [9]. 
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many gluons and light quark-antiquark pairs and 
no particular partonic channel is singled out. In 
the second case we deal with partonic states of 
the type jgg, jggg, jggg, etc., where the photon 
is directly observed in the final state. Since from 
the viewpoint of perturbative QCD a photon is a 
parton as much as gluons and quarks, we are not 
performing an incoherent sum over all the pos- 
sible final states. The final state we are looking 
at is not fully inclusive, and the corresponding 
matrix element M in Eq. (1) is not collinear safe. 

More precisely, the matrix element M(QQ — » 
7 + hadrons) gets contributions from two differ- 
ent components (Fig. 1). Besides the usual direct 
component in which the photon is radiated from 
the heavy quark, there is an additional fragmenta- 
tion component in which the photon is emitted by 
final-state light quarks. The fragmentation com- 
ponent is obviously of leading-twist order (i.e., 
not suppressed by inverse powers of mq in the 
limit mq — » oo) and is not finite in perturbation 
theory. There are indeed singularities from the 
integration region where the photon is collinear 
to the radiating quark. The only consistent way 
of dealing with these singularities is to factorize 
them into non-perturbative fragmentation func- 
tions of the photon D(z, mq). 

Before presenting the complete formalism 
which includes the fragmentation function, let us 
estimate the order of magnitude of the fragmen- 
tation component. The associated collinear sin- 
gularities are hidden in two- (and higher- ) loop 
Feynman diagrams of the type in Fig. 2. Here the 
leading phase-space region is that in which the 
photon-quark angle 9 iq and the quark-antiquark 
angle 9 qq are both small. Correspondingly we ob- 
tain a double-logarithmic singularity of the type 
(lnmg/Qg) 2 , Qo being a collinear cut-off of the 
order of some hadronic scale. Multiplying this 
contribution by the relevant coupling constant 
factor, we obtain aa s (mq) (In niq/Qo) 2 - Since 
lnm^/Qo ~ ^/ a s(mq), we see that the collinear 
fragmentation kills two powers of as thus leading 
to a contribution exactly of the same order as the 
lowest-order direct term, i.e. 0(aa 2 s ). We con- 
clude that the fragmentation component may be 
of the same size as the leading- order direct com- 
ponent. 

This argument also shows that in the domi- 
nant collinear region {0 iq ,9 qq — > 0), the gluon 



which initiates the fragmentation cascade is near 
mass-shell. Therefore the leading fragmentation 
contribution involves D gl (z, mq), the gluon frag- 
mentation function of the photon. Actually, the 
process in Eq. (2) is the only physically relevant 
process in which the gluon fragmentation function 
D gl , and not the quark fragmentation function 
D ql , enters in LO. In fact, in the case of prompt 
photons produced in e + e~ -annihilation only D ql 
is involved in LO, whilst both D ql and D gl enter 
the LO calculation of prompt-photon production 
by hadron collisions. 

As a final comment, we should point out that 
the total width (i.e. the width integrated over the 
energy fraction z = E 1 /mq of the photon) for the 
process (2) is not only collinear unsafe but also 
infrared unsafe. Indeed the fragmentation pro- 
cess has an associated bremsstrahlung spectrum 
of the type dzj z which cannot safely be integrated 
down to z = 0. Only the photon energy spectrum 
(away from z = 0) is computable in (fixed-order) 
perturbation theory, provided collinear singulari- 
ties are factorized into D ql and D gl . 

3. The photon spectrum 

The photon spectrum for the decay of the 
quarkonium state Xqq = J/^,T into 7+ 
hadrons can be written as follows 



dT _ 32 (tt 2 
dz 



9?r 



M 2 dz 1 



(3) 



dT~, 
dz 



■(z, M) = e Q aa s (M) T 7 (z, e Q ;a s {M)) (4) 



-a%(M) 



a=q,q,g 



^ dx z 
— T a (x;a s (M)) D ai {-,M) 



Here M ~ 2 mq is the quarkonium mass and eq 
is the corresponding heavy-quark charge. The 
first and the second term on the r.h.s. of Eq. (4) 
respectively denote the direct and fragmentation 
contribution 2 . Ya {A = J,q,q,g) are hard coef- 
ficient functions computable as power series ex- 
pansions in as'. 



T A (z;as) = T { ° ) (z) + ^T^(z) + 

TT 



(5) 



2 Strictly speaking, the distinction between direct and frag- 
mentation component is factorization-scale dependent. In 
Eq. (4) we have set both the renormalization and factor- 
ization scales equal to M . 
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The direct term to LO is well known [7] and well 
approximated by a linear spectrum (z) ~ 2 z. 
As regards the fragmentation terms f ^ , we know 
from Sect. 2 that only is non- vanishing. Its 
actual computation gives ff ] (z) = 5f^ 0) (z)/12. 

None of the NLO coefficients f ^ is known at 
present. The calculation in Ref. [8] refers to the 
integral dzT~}\z) . As for the shape of the di- 
rect term I\(z) in higher orders, only some model 
calculations are available [11,12]. For these rea- 
sons in the following we limit ourselves to consid- 
ering the photon spectrum in LO. 

4. Leading-order results 

In order to compute the spectrum in Eq. (4) 
we need the fragmentation functions D al (z,M) 
of the photon. Perturbative QCD predicts only 
their M-dependence via evolution equations of 
Altarelli-Parisi type [13]. The solution of these 
equations depends on non-perturbative fragmen- 
tation functions at some input scale Qq. The lat- 
ter are measurable in principle but unknown at 
present. In our analysis we use a reasonable set 
of LO 3 fragmentation functions as provided by 
Owens [13]. This set is obtained by choosing the 
input D ai (z,Qo = A) = 0, A being the QCD 
scale, and is parametrized as follows [13] 

Note that the functions f a (z) are scale invariant, 
so that the full M-dependence is due to the factor 
in front of them on the r.h.s. of Eq. (6). The 
actual form of f a (z) given in [13] shows that the 
gluon (quark) fragmentation dominates at small 
(large) values of z. 

The LO photon spectra obtained by inserting 
Eq. (6) into Eq. (4) are presented in Fig. 3. Note 
that, because of the dependence D gl oc aj as, 
the ratio between direct and fragmentation con- 
tribution is independent of as(M) (the value of 
as(M) only affects the overall normalization in 
Fig. 3). Note also that this ratio is proportional 
to e'q (see Eq. (4)), so that the relative effect of 
the fragmentation is larger for T than for J/tp. 



3 NLO fragmentation functions are also available [14], but 
they are not relevant in our LO analysis. 



The modification of the LO direct spectrum 
due to the fragmentation component reflects the 
z-behaviour of the gluon fragmentation function. 
The latter is very soft and, hence, the fragmen- 
tation contribution is dominant at small z but 
quite small at large z. On the other hand, ac- 
curate experimental data [15,16] are available at 
present only in the large-z region. It follows that 
the LO fragmentation term considered here does 
not have a big impact on the QCD studies per- 
formed so far (typically, only the data points for 
z^0.5 have been used in these analyses). 

The same conclusion, however, does not nec- 
essarily apply beyond LO. In fact, although the 
NLO fragmentation term is still unknown, also 
the quark fragmentation function of the photon 
is involved in this order. According to the Owens 
parametrization [13], D ql is likely to be two or- 
ders of magnitude larger than D gl at z^0.6. 
Therefore the NLO fragmentation component is 
not necessarily small in the large-z region. More- 
over, the NLO direct contribution to dT/dz is also 
unknown and, more importantly, one cannot con- 
centrate on very large values of z in the attempt 
of suppressing the effect of fragmentation. In the 
'extreme'-z region (say, 1 — z^(lGeV)/M, i.e., 
z^0.7 for J/tp and z^0.9 for T), hadronization 
effects are much larger than contributions com- 
putable in perturbation theory. 

As a result of this discussion, we regard the 
present tests of perturbative QCD (in NLO) from 
the process (2) as quite uncertain. 

5. Conclusions 

We have shown that the perturbative QCD ap- 
proach to the radiative decay of quarkonia re- 
quires the introduction of a fragmentation compo- 
nent (surprisingly never considered so far) besides 
the usual direct contribution. Such fragmenta- 
tion component enters to leading-twist order and 
to LO in as- The perturbative QCD predictions 
for the photon spectrum thus depend on as and 
the fragmentation functions D al of the photon. 
In principle, measuring the photon spectrum at a 
single scale ( J/tp or T decay), one cannot extract 
a value for as- Measurements from both J/tp and 
T decays, in contrast, allow a combined determi- 
nation of as and D ai . In practice, the fragmenta- 
tion contribution is small at large z, so that data 
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in this region can be used for a^-determinations. 

On the other hand, fragmentation dominates 
at small values of z. In this region it is experi- 
mentally very difficult to disentangle the prompt- 
photon component from the large background due 
to 7r and rj decays. Nevertheless, further exper- 
imental and phenomenological investigations of 
this issue can be extremely interesting. Data at 
small z could indeed be used for measuring (or 
setting bounds on) the gluon fragmentation func- 
tion D gl of the photon. The perturbative QCD 
predictions for prompt-photon cross sections in 
hadron collisions substantially depend on D gl , 
and at present we do not have any experimental 
information on it. 

As for a^-determinations from quarkonium de- 
cay, we would like to add a last comment. The 
ratio between the hadronic and leptonic widths 
is known up to NLO in perturbative QCD but, 
due to the different topologies of the hadronic- 
and leptonic-decay mechanisms, it is likely to be 
affected by large relativistic corrections [4,5,17]. 
These corrections are expected to be smaller for 
the ratio between the radiative and hadronic 
widths, but, as discussed in this paper, a com- 
plete NLO calculation for the photon spectrum 
is still missing. For these reasons we recommend 
(see, for instance, Refs. [5,18]) a careful and 'gen- 
erous' estimate of the theoretical uncertainties on 
the values of as extracted from quarkonium de- 
cays. 
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Fig. 1. (a) Direct and (b) fragmentation con- 
tributions to the matrix element M(QQ — > 7 + 
hadrons). 

Fig. 2. The lowest-order Feynman diagram 
leading to photon fragmentation. 

Fig. 3. LO results (in units Rq/M = 1: see 
Eqs. (3) and (4) ) for the prompt-photon spec- 
trum in (a) J ftp and (b) T decays, with as = 0.25 
and as = 0.20 respectively. 
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